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THERMAL NOISE IN SPIN-PHONON SYSTEMS 


Harold A. Sabbagh 
Rose Polytechnic Institute 

ABSTRACT 

Starting with the combined lattice and spin- lattice Hamiltonian densities, 
field equations for transverse elastic waves are derived. An electrical trans- 
mission line analog is established which employs an "impedance" function 
(matrix) calculated using the density- matrix formalism for the interaction of 
a paramagnetic spin system with acoustic vibrations. 

Associated with the impedance function is a continuous distribution of 
thermal noise "voltage" sources whose statistical characteristics are de- 
termined by a form of the (Nyquist) "fluctuation-dissipation" theorem. The 
significance of positive and negative temperatures on the spectral-density of 
the noise sources is pointed out. 

Application of the theory is made in determining the spatial correlation 
and spectral densities of the resultant acoustic noise fields in infinite and 
bounded systems. 
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THERMAL NOISE IN SPIN-PHONON SYSTEMS 


INTRODUCTION 

With the advent of acoustical devices in engineering systems/ 1 ) it is de- 
sirable to investigate noise properties of phonon systems. A particularly in- 
teresting study is that of paramagnetic ions doped into an elastic medium, for 
this can serve as a model for a spin-phonon maser, adjustable acoustical delay- 
line, etc. 

The problem of thermal noise in electromagnetic systems has been in- 
vestigated (Haus, 1961; Vanwormhaudt and Haus, 1962; Landau and Lifshitz, 
1960), and has been applied to electromagnetic masers (Siegman, 1964). It is 
appropriate to adapt some of the ideas presented in those earlier studies to the 
problem at hand, namely thermal noise fluctuations in an elastic system. 

Our approach to the problem is to determine the lattice and spin-lattice 
Hamiltonian densities and then derive the canonical field equations appropriate 
to the model. We then deduce the spin-phonon susceptibility matrix which re- 
lates the strain-vector to a spin-induced stress vector and from this determine 
an electrical transmission-line analogy. 

At this point contact can be made with the earlier investigations mentioned 
above in determining the spectral densities of certain fictitious noise "voltage" 


(^See for example Proc. IEEE (special Issue on Ultrasonics), Vol. 53, No. 10, October 1965. 
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sources. These sources act as inhomogeneous terms in the transmission-line 
equations, which are solved for either unbounded or bounded configurations 
using Green’s functions. Once the "voltage" and "current" on the transmission- 
line are determined.it is quite straightforward to calculate their spectral-densities. 

HAMILTONIAN DENSITY FOR A SPIN-PHONON SYSTEM 

Hamiltonian Density of an Elastic Medium 

The Hamiltonian density of an elastic medium is an expression of energy per 
unit volume given in terms of momenta, displacements, derivatives of displace- 
ments (strains), etc. The total energy (kinetic plus potential) residing within a 
certain region of space is given by the integral over that region of the 
Hamiltonian density. 

The lattice (elastic) Hamiltonian density for an isotropic medium is given 
by (Morse and Feshbach, 1953): 



where p - mass density, and p. and \ are the two elastic constants specifying the 
stress-strain relationship of an elastic medium. u x , u , u z are the components 
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of the vector displacement of the medium at x, y, z, while p x - p(d u^dt), 

P y = p(dy/dt) t p z = p(du z /dt) are the components of the momentum density 
at x, y and z. 

The momentum terms in (1) are the kinetic energy density while the re- 
maining terms are the potential energy density arising from the stress-strain 
interaction. •* 

Hamiltonian of a Paramagnetic Spin System (Spin Hamiltonian) 

The spin Hamiltonian is a quantum mechanical operator function of spin 
operators (Bleaney and Stevens, 1953; Bowers and Owen, 1955). Its eigenvalues 
are the allowed energies of the (unperturbed) paramagnetic ion, and its eigen- 
vectors are the states corresponding to the allowed energies. 

The energy levels are important to the determination of operating (resonant) 
frequencies of paramagnetic devices, such as masers, while the eigenvectors 
are used to determine the matrix elements, or transition probabilities, which 
indicate the "strength" of interactions between the elastic system and the spin 
system. 

Interaction Hamiltonian Density (Spin-Lattice Hamiltonian) 

In order to describe the interaction between the elastic and spin systems, 
one needs a Hamiltonian density which involves elastic variables (displacement, 
strain, etc.) and spin operators (Tucker, 1965). There are three important con- 
ditions that such a Hamiltonian must satisfy (Van Kranendank and Lee, 1966): 

/ 
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1. Hermiticity : The spin-lattice Hamiltonian density must be Hermitian. 
This guarantees that the interaction energies (the eigenvalues of the 
Hamiltonian) will be real. In order that the Hamiltonian be Hermitian 
it must be a real function of the elastic variables, d.c. magnetic field, 
and spin operators. All combinations of spin operators appearing in 
the Hamiltonian must, themselves, be Hermitian. 

2. Time-reversal invariance : Because the spin-lattice coupling involves 
electromagnetic interactions, due to the presence of a spin-vector and 
an external magnetic field, the spin-lattice Hamiltonian must be in- 
variant to time- reversal. This means that upon changing the signs of 
the external d.c. magnetic field and spin- vector, the Hamiltonian becomes 
complex-conjugated. 

3. Point symmetry : Finally, the spin-lattice Hamiltonian must be in- 
variant under the group of operations corresponding to the point sym- 
metry of the paramagnetic ion's site in the host crystal. Thus, if we 
assume that the ion occupies a site of octahedral symmetry within the 
crystal, the interaction Hamiltonian density must be invariant under 
the operations of the octahedral group. 

These conditions imply that the continuum (long-wavelength) limit there 
will be four possible forms for the open-lattice Hamiltonian density, 3^: 


Kl x = ^a x (h 0 • s)div (u) 
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= /^A 2 [(H 0 S + SH 0 ) : e- J(H 0 -S)div(u)] 

= A'N A 3 ( H 0 x s ) • curl (u) 

3 

= /jNA 4 [ss : 1 - J S(S+ l)div(u)] . (2) 

Aj - A 4 are empirically determined constants, N = no. of ions per unit 
volume, H 0 is the d.c. magnetic field, and S is the ion's spin-vector operator. 
Also e ab = 1/2 (<9u b /<9x a + <9u a /dx b ) is the strain-tensor and 

SS • T = 2Z S * S b e b. • 

a,b 

The first three forms in (2) are linear in each of the variables H 0 and S, 
whereas the fourth is independent of H 0 but depends quadratically on S. Con- 
trast these possibilities with the single Hamiltonian, g/3H • S, of electromagnetic- 
spin interactions. 

The total Hamiltonian density for the combined spin-elastic system is given 
by the sum of the three Hamiltonians mentioned so far: 

U = U, + u s + U A , ( 3 ) 

where U s is the spin- Hamiltonian. 

As (3) stands it is a mixture of a classical Hamiltonian, and quantum- 
mechanical Hamiltonians, and # s ,£, involving spin operators. Eventually 
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we will take quantum average values of these operators and proceed 
classically. 

CANONICAL FIELD EQUATIONS 

One important reason for expressing the mathematical model of a physical 
system in terms of a Hamiltonian density is that one may derive the differential 
equations for a variety of different subsystems of the given system in a canoni- 
cal form. 

Because the spin subsystem in our problem requires a quantum- mechanical 

Hamiltonian formulation, it is natural to develop the elastic subsystem in a 

Hamiltonian manner, also. As previously mentioned, the latter Hamiltonian is 

classical, i.e., it does not involve non-commuting operators. 

From here on we consider transverse vibrations propagating in the 

z -direction, the assumed direction of the d.c. field H 0 . Thus, in (1) and (3) 

we set p z , d/dx, d/dy, u z , H 0 and H 0 equal to zero. Note also that terms in- 

* y 

volving div (u) vanish. The resultant combination, H', of the lattice plus spin- 
lattice Hamiltonians becomes, after taking quantum- average values 

H' = ip(P**P y 2 )*? 

<?u <9u 

■ + -J7 + -d • (4) 
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where 


M x = <(a 2 + a 3 )H(,s x + a 4 (s z s i + s x s,)) y 

M y - <( A 2 + A 3 ) H oS y +A 4 (S,S y +S s S,))rtJ , (5) 

and the brackets denote quantum- averaging. 

Hamilton’s canonical equations are [Morse and Feshbach, 1953] . 


<?H' 

a aH' a 

aH' 

a 

aH' 

aH' 

u x - <9p x * 

Px Ax (du \ + Ay 

Aw) 

(du\ - dz 

Aw) 


au 

X 

aH' 

a aH' a 

a h' 

a 

aH' 

an' 

u y " <3p y ; 

p y ax 7au \ + dy 

Adi) 

/au \ dz 

a (ar) 

e 

au 

y 

air 

a aH' a 

aH' 

a 

aH' 

aH' 

u z " a Pz s 

p z (9x /Au \ + <?y 

H"aw 

(du \ + dz 

d \w) 

<(S 

Au 

z 


where a dot denotes time-differentiation. 


Upon substituting (4) into (6), the following system is obtained 


u 


“ P x/P 


P 


d 2 u 

X 

dz 2 


dM x 

+ dz 


U y - Py/P i 


P 


d 2 U 

y 

<9z 2 


+ 


-AM 

y 

Az 


( 6 ) 


(7) 
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All terms in (7) are classical; quantum-mechanical spin operators have 
been replaced by their average values. Equation (7) shows that the source of 
elastic vibrations is the gradient of the paramagnetic stress functions M x and 
M y . This agrees with the statements of Jacobsen (1960) and Seavey (1965) re- 
garding piezoelectric and magnetic production of acoustical vibrations, 
respectively. 

M x and M y are not given a priori . They depend on Ju x /Jz and du y /dz (actu, 
the randomly fluctuating parts of U and M„ can be considered to be given func- 
tions, independent of du x /dz and du y /dz . We show this in another section). 
After determining M x andM y as functions of <9u x /<9 z and <9u y /<?z, we substitute 
the results back into (7) to obtain wave equations for u and u . The relation- 
ship between M x , M y and du x /dz , du y /dz involves the spin-phonon susceptibilit 
matrix. 

SUSCEPTIBILITY MATRIX 

Because the spin system is an atomic system we must resort to quantum 
mechanical principles for calculating average values of the pertinent spin 
operators. A useful tool for this purpose is the density matrix (Sabbagh, 196 
Slichter, 1963). 

Every quantum- mechanical operator may be given a matrix representat: 
by choosing a complete, ortho- normal basis system of vectors. The basis tl 
we will use will be the eigenvectors of the spin-Hamiltonian. 
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The average value of operator. A, then is given by 

(A) = Tr[^..][A jn ] , (8) 

where Tr denotes the trace (sum of diagonal elements) of a matrix, p i j is the 
ij — element of the density matrix, and A - n is the jn— element of the matrix 
representation of the operator A. If there are p energy levels in the paramag- 
netic spin system, all matrices will bePxP, 

To a first (linear) approximation the diagonal elements, p mm , m = 1 . . . , p 
of the density matrix remain constant at their thermal- equilibrium values, pJJ~> . 
The off-diagonal elements, to the same order of approximation, satisfy the 
system of linear equations 

~ST = (i W mn - rM P mn + TT ~ PnP) U U m - ( 9 ) 

where m / n and m, n = 1, . . . , p. a J mn = E n - E^/R is the frequency correspond' 

ing to the difference between the n— and m— unperturbed energy levels. The 
r*s are time-constants related to the return of the spin system to thermal 
equilibrium after a transient perturbation. Finally, U p is the mn— matrix 

mn 

element of the spin-lattice Hamiltonian operator. All matrix elements are 
between the unperturbed eigenstates of the spin-Hamiltonian. 
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In the sinusoidal steady-state & s/ £ , for a single ion, becomes (see (5)) 

mn 



/ au x 

au 

\ 



) = ls x ' A* "37 + 

ft dz 

e h "t 


mn 

\ mn 

y mn 

/ 



( < 


aU y*\ 



+ V x mn ^ dz 

\ mn 

+ s; v 

y mn 

az Je-- 

(10) 


The asterisk denotes complex conjugate and the field amplitudes U x , U y are 
independent of time. In (10) we have made the definitions 


S x ' 

= (a 2 + A 3 )h,,S, * A 4 (8,S,+S,S.) 


s ; 

= (A 2 + A 3 )H 0 S y + A 4 (S Z S y + S y S z ) . 

(11) 


The solution of (9), given (10), is easily obtained because of the sinusoidal, 
steady- state condition (it is equally easily obtained as a convolution integral, 
as well): 


j( (T) _ <t>\ 

J \'mm r nn / 


v- 


p mn (*) 


au 


(9U 


S' 


■U dz + S 4m dz 

D- 


L 


cot 


(s. 

\ x ^ 


au 


au 


U + S 4 m 

D+ 


-) .1 

— e" J"t 


( 12 ) 
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where 





1 

Hm 


D+ - -}(^ m +co) + — • (13) 

Comparing (5) and (8) we see that 

M * = 

M, = . (14) 

After carrying out the operations implied in (14) (using (12), of course), 
dropping any time- independent terms and rearranging, we get 

5U X . «3U y 

M X a xx M dz eJa>t + a xy M Tz eJ£Ut + com P lex conjugate 


5U x . flU y 

M y <9z eJ&,t + a yy M TT eJ<at + complete conjugate , (15) 

In (15) 

a xx == IT 2Z ( n i <T) ~ n m (T) ) 

l<m 



a 


xy 




(16) 
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a 


yx 


Km 


S y jp S xV 

*mii X/ Vm 

EF 


o * o • 

■ X U Hm 

' “ D+* ~ 


yy 


TT ( n l <T)_n « (T> ) 

! <m 


q # q / q * O # 


D- 


D+ 


( 16 ) 


= thermal equilibrium no. density of ions in state 1, etc. 

The two-by-two matrix with entries a. , a , a. and a is the spin-phonon 

xx x y yx yy 

susceptibility matrix. 

Consider a xx , for it is typical of the other three functions. It consists of a 
summation over pairs of energy levels (with l<m) of terms whose real and 
imaginary parts are of the form 


Real part 


Imaginary part 


|S4I 2 


— 


Hm Km " 

» 2 ) 

Km "“ 2 ) 2 + 

4co 2 
_ 2 




2 fjL 

K 






c 2 ) 2 


Aw 2 

+ 7 J 


(17) 


We see quite clearly the role played by the matrix .elements S/ m in determining 
the interaction strength between spins and phonons . 

The functions shown in (17) are remiiiiscent of the real and imaginary parts 
of immitance functions of RLC circuits. w£ ro is the resonant frequency of the 
resonant circuit corresponding to the ^m th pair of energy levels and 2/v^ m is the 
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half-power bandwidth of this same "circuit.” We conclude, therefore, that the 
spin-phonon susceptibility functions are analogous to a series combination of 
RLC resonant circuits. 

ELECTRIC TRANSMISSION-LINE ANALOGY 

To develop the electric transmission- line analogy we start with (7), which 
may be combined to give 

* 2 “, u, 


a ! u y _ 3*u y £M y 


We establish the analogies 


<9u x <9u x 

<9t ^ ^ dz ^ V x 


(18) 


du <?u 

TF ~ i y 5 7 " V y • ( 19 > 

Note that when fidu/dz (= stress) is multiplied by <?u/<9t (velocity), power/ 
area is obtained. If we imagine one-dimensional wave propagation through a 
unit cross-section, then stress x velocity x ^m 2 yields power which is v x i. 
Hence, the aptiiess of the analogy. 
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In terms of the "electrical” variables, (18) becomes 


^X 


1 ^x 

di 

X 

p dt 

dz + dz ’ 

fJi dt 

= ST 

di 

y 

<?v y dM y 

1 dv 
± y 

di 

y 

Pit 

dz + dz ' 

ji dt 

“ dt ' 


( 20 ) 


Upon taking the Fourier transform of (20) and using (15) we get the vector- 


matrix equations 



V 


ri+a 

XX 

a xy 


~dV x /d z 

j up 

- 1 *- 


a 

L. yx 

1+ %y- 


_dV y /d z _ 


“v x “ 

1 

j 

"dl x /d“ 



yup 

_ V y_ 


_d! y /d z _ 




To simplify matters, we assume a wave linearity polarized in the x-direction. 
Therefore, after setting V y = I y = 0 and dropping subscripts, (21) becomes 


; <*P T 

j 1 + a * 


dV 

d 


co 

j 7T V 


DI 


( 22 ) 
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These equations are obtained from the incremental transmission-line con- 
figuration shown in Figure 1(a). Observe the series impedance 

j^Az _ i^ Az ( i + a R -j a l) f 
1 + a 1 1 + a| 2 

where a R is the real-part and a x the imaginary-part of a, and shunt admittance 

j Ct)/ (A Az. 

When we recall the summations appearing in (16), we may redraw Figure 
1(a) using a series combination of impedances as in Figure 1(b). Each im- 
pedance is associated with an energy-level pair as in (17). 

Now that we have established an electrical analogue, it is rather straight- 
forward to develop the thermal noise source associated with the dissipative 
mechanism of the transmission-line. 



Figure 1(a). Incremental Section of Analogous Transmission- 
line Yielding Equations (22) 


15 





Figure 1(b). Illustrating Expansion of Impedance into a Series Combination in Accordance 
With (16). Each Block Corresponds to an Energy-level Pair. 

THERMAL NOISE SOURCE 

An expression for power-dissipation in terms of a may be obtained from 
(22) by multiplying the top equation by I*, conjugating the bottom, multiplying 
the result by V, then adding the two equations and taking one-half of the real-part: 


dP 

Hi 


2 cop ct-j 
2 I 1 + ct | 2 


III 2 



i= 1 


(23) 


The second equality in (23) shows that each energy-level pair contributes 
to the overall change in power along the line. If a x is negative for the i— 
energy-level pair, then that pair contributes gain (maser action at the correspond- 
ing resonant frequency). Otherwise we get dissipation' (or absorption) of power. 

Now we shall show that a x is associated with the spectral density of noise 
sources arising from the series impedances. The generalized Nyquist theorem 
(Davenport and Root, 1958; Twiss, 1955) asserts that the thermal noise associated 
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with a lossy network at temperature T is equivalent to fictitious noise sources 
connected to the terminals of the network. The spectral densities of the sources 
are related to the dissipative properties of the network. This result, general 
as it is for electrical networks, is a special case of the "fluctuation-dissipation" 
theorem satisfied by all physical systems in thermal equilibrium (Landau and 
Lifshitz, 1958; de Groot and Mazur, 1962). 

By the spectral density of a random process, x(t), we mean 

S x («) = limE[|X T ( w )| 2 ], (24) 

where 

1 f T 

X t<» = Tff J x ( t ) e~ J a>t dt , (25) 

and E[ 3 implies the expected value of whatever is within the braces. Through- 
out this report when we speak of the Fourier transform of a function we mean (25). 

The significance of the spectral density lies in the fact that the expected 
mean-square value of x(t) is given by 





(a)) d co . 


(26) 


By the "one-sided" spectral density we mean 2 S x (oj) with the integration in (26) 
extending only over positive frequencies. 
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Before determining the noise sources associated with each block in Figure 
1(b) we make the following preparatory remarks. The noise properties depend 
on the "temperature" of each block. If the spin system were in thermal 
equilibrium with the lattice, each energy level would be populated according 
to the Boltzmann distribution with a single temperature T. Then the ratio of the 
populations of any pair of energy levels would be simply 

n E. - Ej 

nT = — kT- (27) 

where k is Boltzmann’s constant and E. and E. are the energies of the i — and 
j — levels . 

In the presence of an external source (such as a pump in a phonon maser), 
however, the distribution of populations will be altered radically. For example, 
one pair of energy levels may have an inverted population (more ions in the 
higher than lower energy level), and another may have equal populations etc. 

If we assume that we may always define an equilibrium temperature between 
pairs of energy levels so that (27) is satisfied, then there will be as many 
equilibrium temperatures associated with the spin system as there are pairs 
of energy levels. In the examples above the inverted populations correspond 
to a negative temperature, while the equally populated pair of energy levels 
corresponds to infinite temperature. 

An important result is that there is a correspondence between a x and 1\, 
the temperature of the i — pair of energy levels (or the i— spin subsystem). 
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When T. > 0, (17) and (27) show that a z > 0; T* > 0 corresponds to a r < 0, and, 
finally, when T £ = then a.j * 0, These results hold for positive frequencies. 
Because a : is an odd function of oj, the reverse of each holds for negative 
frequencies. From (23), however, we note that the significant function is coo. l , 
which is positive, negative, or zero, corresponding to the above temperature 
ranges, at all frequencies. 

In order to determine the spectral densities of the fictitious noise sources 
associated with each impedance block, we remove the m— box in Figure 1(b) 
and connect a series noise voltage source to account for the thermal noise fluc- 
tuations at temperature T m (see Figure 2). To this combination we connect a re- 
sistor R m Az at temperature T m (Figure 3(a)). 

The spectral density of noise power delivered to R ra Az is easily calculated 
to be 


Sr 


1 

2 


S tn O) R m AZ 

(R. + zJ(«. + C) ' 


(28) 


where S m (co) is the spectral density of e m . 

The resistor at temperature T m is, itself, a noise source and may be con- 
sidered to be a voltage source e R Az in series with R ra Az (see Figure 3(b)). The 
spectral density of the noise power delivered by the resistor to Z M Az is 


S 


z 


! S. R MRe(z„)Az 

1 ( R . + Z J( E „, + z .*) ' 


(29) 
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At thermal equilibrium (28) and (29) must be equal, which means that 


S e («u) cu pa. 


S m<» - R~ 


1 + a 


( 30 ) 








Figure 3(b). Circuit for Calculating Resistance Noise Power 
at Temperature T m Delivered to Z m Az. 

The one-sided spectral-density of e R Az is (Davenport and Root, (1958)) 
4kT m R m Az, which implies that S e (w) = 4kT in R m /Az. Thus, (30) becomes 

4kT w P a i 

s„o> = -sr ' ]777ji • < 3l > 

This is the spectral density of a source at temperature T m lying within the 
interval (z, z + Az). Such a source is statistically independent of all sources 
lying within any other intervals, no matter how fine we make such intervals. 
Thus, taking the limit Az - 0 in (31) we obtain for the expected value of 

e mOl) * e m*( Z 2> : 

cop a j 

E [ e m( Z l) * C ( Z 2 )] = 4kT m | { + J 2 8 ( Z 1 " Z 2 ) > ( 32 ) 

where 8 (z x - z 2 ) is the Dirac delta function. 
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kT m is the low-frequency (or high-temperature) limit of the quantum- 


mechanical expression 



which means that our final result for the expected value of e m (z x ) • e m * (z 2 ) is 



This is the important relation between a x and the noise source at temper- 
ature T m . Note that the frequency dependent coefficient of the delta function in 
(33) is always positive, regardless of the sign of T m , because, as was shown 
above, T and a. take the same sign. Another important feature is that (33) 
has a peak at the resonant frequency of the m th energy-level pair and then falls 
off to essentially zero outside the 2/r bandwidth centered at the resonant frequency. 

If the i th and j th energy levels together comprise the m th energy- level 
pair and if the no. density of ions occupying this pair of energy levels inN m then 


n i + n j 


N 


m 



( 34 ) 
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imply that 



where K&> = E. - E.. . 

m j 1 

Upon comparing the imaginary part in (17) with (33) and (35), we see that 
(33) is proportional to 



4K« 



As T m approaches infinity this function approaches 2N m E«; it approaches zero 
as T m -0 + , and, finally, it approaches 4N m Hoj as T m - 0_ (complete population 
inversion). Each of these results is consistent with the fact that thermal noise 
of this type originates with ions occupying the higher energy level of a given 
pair. Such spontaneous emission noise is discussed with reference to electro- 
magnetic masers in Siegman (1964). 

Each generator associated with an impedance block in Figure 1(b) (see Fig- 
ure 2) satisfies (33) with the appropriate temperature and a T . In addition such 
generators are statistically independent of each other, even though they lie 
within the same interval. Because of the series connection of impedance blocks 
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and noise sources and because of the independence of these sources, it follows 
that the total noise voltage source, eAz, is just the sum of each of the individual 
noise sources within the small interval and we have 

E [ e ( z i) • e *( z 2 )] = 2Z E [ e m( z i)* e * ( z 2 )] > (36) 

m= 1 

where P, again, stands for the total number of impedance blocks (which is equal 
to the total number of energy- level pairs in the paramagnetic spin system). 

If we insert eAz in series with the series impedance in Figure 1(a) we can 
derive an inhomogeneous system of equations satisfied by the fluctuating field. 
Kirchoff's voltage and current laws give, respectively, 

V(z + Az) = I(z+Az) ^ - + eAz + V(z) 


I(z + Az) - I(z) + j Az V(z) . 


(37) 


After rearrangement, division by Az and passage to the limit Az - 0, we get 


dV 

Hi 


I 


i up 

1 + a 


e(z) 


dl 

Hi 



V(z) . 


( 38 ) 


Equation (38) differs from (22) only in the explicit inclusion of the com- 
posite noise source voltage-gradient e(z). 
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Finally, let us say a few words about the spectral density of the random 
sources associated with the complete vector system (21). We will consider the 
special case in which a xx = a yy = a t , a xy = -a yx = a 2 , because these con- 
ditions yield circularly polarized acoustic waves and have been the subject of 
some study (Muller and Tucker, 1967). 

The space rate- of- change of time-average power along the transmission- 
line is given by 



where A is the matrix appearing in (21), A -1 is its inverse, and the superscript 
dagger denotes conjugate-transpose. 

By the assumed direction of current flow (right-to-left), the right-hand side 
of (39) must be non-negative for any I x , I y if the transmission- line is dissipative 
(absorptive). A straight-forward, though tedious, calculation shows that 
j &P [a->-(a- 1 ) + J is indeed a positive, semi- definite matrix for a normal 
(non-inverted) population. 

Because we have a vector system we must associate two random functions 
e x (z) , e y (z) with the dissipation matrix j cap [a -1 - (a _ 1 ) + ] just as we as- 
sociated a single source with the impedance joip/ 1 + a. Corresponding to (33) 
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we have the spectral density matrix 


e x( Z l) e x*( Z 2 ) e x( Z l) e y*( Z 2 ) 

sw*;w 


4Ew 


exp 


m- 


[*-*-(*-*)*] , 

2 ^ ( z 1 — z 2 ) > (40) 


* J up 


which is positive semi-definite. 

The inhomogeneous equations corresponding to (21) are 


_d 

dz 


■v; 

V 

- j cap A 1 

V 

I 

+ 

e x 

e 

y 


y 


y 


_d_ 

dz 



(41) 


CALCULATION OF OUTPUT NOISE SPECTRAL DENSITIES 

As applications of the preceding developments we shall calculate the "voltage” 
and "current" (actually, strain and velocity) spectral densities for various con- 
figurations. The first configuration will be simply an infinite slab of unit cross- 
section, and we shall determine the spatial correlation function as well as 
spectral density. The other configurations will both be finite systems with a 
signal source at one end and the other end either free or terminated in the ap- 
proximate characteristic impedance of the acoustic medium. In either of the 
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latter two cases we will calculate the output spectral densities. The solutions, 
of course, are determined as boundary value problems. 


Correlation and Spectral Density of Noise in an Infinite Region 
We start with (38) expressed in vector-matrix form 


dV 

3z 


MV + e 


(42) 


where 


V 




A 

cop 



-V' 

, M - 

u 

i 1 + a 

j ^ 

"e" 

.1- 


CO 

J 7 

0 


.0. 


(43) 


The solution of (42) is written in terms of a 2 x 2 Green's matrix, G(z - z ' ): 


V(z) 



G( z - z ' ) e( z ' ) dz ' . 


(44) 


where g(z - z ' ) satisfies the matrix differential equation 


dG 

di ~ MG 


8(z - z')U . 


(45) 


U is the 2x2 identity matrix. 



In component form (44) is 


V(z) - g 1]t (z- z') e(z') dz' 


(46) 


1(0 


■r 


’21 


(z - z')e(z') dz' 


g xl and g 21 being elements of G. 

After some manipulations and use of (33) and (36) we obtain for the spatial 
correlation functions of V( z ) and I(z): 


MO s E 


[ V ( z i) V * (z 2 )] “ (J^ gn( r + 0 gfi(Od^ 


p pa)a, 

4K a> T m 


L 


m- 1 exp 


te )-. 1 - 1 


(r) ■ E [l(zj) I* (z 2 )] - g 21 (rH) g 2 *(Od[J 


P 



m=l exp 


4Hoi 



1 


pOOCLj 
1 m 


|l+a | 2 


. (47) 
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where r = z x ~z 2 . These results show that the random processes V( z) , I(z) 
are stationary in space in the sense that their spatial correlations depend only 
on the difference r = z x - z 2 and not an z p z 2 , individually. 

From here on we will assume that all energy- level pairs are at the same 
temperature so that the summation appearing in (47) is replaced by the single 
term 



pa jo-j. 

1 + a| 2 


SO) . 


Because the V and I processes are stationary it is convenient to introduce 
the wave-number representations, H v (k), Hj (k), in order to calculate the 
correlation functions: 

H v (k) 

H r (k) 

Upon substituting (47) (which are really convolution integrals in space) into 
(48) we obtain 

H v (k) = S( ffl )O n (k)G*(k) 

Hj (k) = S.(&>) G 21 (k) G* x (k) (49) 


= f 

j -00 


) e" jkr dr 


■r 


R (r) e _jkr dr . 


(48) 
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where G lx (k) and G 21 (k) are spatial Fourier transforms, respectively, of 
gji (z), g 21 (z) (the wave-number representation, (48), is, of course, the 
spatial Fourier transform). 

The calculation of G n (k), G 21 (k) proceeds most simply by taking the 
(spatial) Fourier transform of (45): 

(jkU-M)G(k) - e" jkz 'u (50) 

G(k) is then given by 

G(k) = e -jkz ' (jkU-M) -1 , (51) 

which, after substituting the matrix M, from (43), and identifying the 11 and 21 
elements of the result, yields 


G u (k) 


jke~ jkz ' 


-k 2 + 2 



G 2 l( k ) = 




( 52 ) 
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Upon substituting (52) into (49) and then taking the inverse transform of (48), 


we get 


R v ( r > = 



k 2 e j kr dk 

“ 2 (i) TTa - [“ 2 (£) TTai- k \ 



We use residue theory (Churchill, 1960) to evaluate these integrals. The 
poles of the integrands are located at a + jb, a - jb, -a - .jb, -a + jb in the com- 
plex k- plane, where 



Note that the poles are symmetrically distributed about the real and imaginary 
axes. 

We close the contour with a large semi-circle lying in the upper half-plane 
for r > o and lying in the lower half- plane for r <0. This assures that contributions 
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to the overall integral are negligible as the radii of the semi-circles go to 
infinity. The calculation of the residues is quite straight-forward with the re- 
sult that 


_ S(oi) . , , I , 

R v ( r ) - e | r| (a cos ar -b sin a | r| ) 


MO 


S<«) 

4ab 


O’ 


e“ b I r 


t 


a cos ar + b sin a 


a 2 + b 2 


-) 


(55) 


The correlation functions are even in r. We introduce the terms correlation 
length and correlation wavelength for 1/b and 277 /a, respectively. Due to the 
presence of an imaginary part in a, we see that the farther apart in space two 
noise samples are measured (in the frequency domain) the less correlated they 
are. 

Upon letting r = o in (55) we obtain the spectral densities S v (oj), Sj (oi) of 
the V and I processes, respectively. Thus, substitution of (54) into (55) gives us 


s v O) - 


Sj<*0 = 



(56) 


32 




These results admit the following interesting interpretation. From Figure 
1(a) the "inductance" per-unit-length is seen to be p / 1 + a and the "capacitance" 
per-unit- length is -l/p. Hence, the magnitude of the wave impedance is 



(57) 


Therefore, upon dividing S v (w) by (57) or multiplying Sj (a>) by the same 
function we get the spectral-densities of the power associated with the "voltage" 
and "current" processes, respectively. The power spectral- densities are equal 
and their sum is 


S p («) •= 



(58) 


If the spins are brought into thermal equilibrium with the acoustic medium 
at temperature T, and then if the medium is made nondispersive by setting 
a(co) = 0, we get the phonon thermal equilibrium result 


S p <»> 


4Ho) 

/fTco\ 

exp \ W ) - 1 
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Our results are general in that they contain this result as a special case. The 
corresponding result for electromagnetic radiation in 3-dimensional space is 
given in Landau and Lifshitz (1960). 

Signal Source and Free Termination 

Consider a signal source connected at z = 0 and the end z = L left free. If 
we take the actual (acoustic) source to be an imposed strain then the electrical 
analogue is a voltage source. The free surface at z = L appears electrically to 
be a short-circuit because a free surface cannot support a shear-stress. Since 
stress is given by 

<9u <?M du 

P dz + dz “ K 1 + a ) Jl ~ ( 1 + a) V , 

we see that V(L) = 0, hence, the short-circuit termination. 

Our concern is with the calculation of the spectral-density of the current 
(velocity) at z = L due to the signal and noise. We take the signal and noise to 
be uncorrelated. We start with (38) and the boundary conditions 

V(O) = E s (a;) 

V(L) = 0 , (59) 

where E s (a>) is the Fourier transform of the "voltage" signal source connected 
at z = 0. 
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The output current will be a superposition of the currents arising from E s 
and the noise source e(z), each acting alone. We label the voltage and current 
arising fromE s , V s and I s , while that resulting from e(z), V n and I n . V s and I s 
satisfy (38) with e(z) set equal to zero and the boundary conditions (59). V n , I n 
satisfy (38) and the boundary conditions V n (0) = V n (L) = 0. Then V = V n + V s 
and I = I n + I s satisfy (38) and (59). 

It is a simple matter to show that V s and I s are given by 


where 


V s (z) 


E s (“) 


sin k (L - z) 
sin k L 


M z > 


. f cos k (L - z) 
\Z 0 y sinkL 


(60) 


k 


CO 


( m ( 1 + <*)) 


1/2 


and 


Zo 



1/2 


We next postulate a noise solution of the form 


■f 

Jo 


V n (z) - g y (z, z') e(z')dz' 


In(*> 


-^0 


z') e(z') dz' 


(61) 
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with the voltage and current Green's functions satisfying 


dz J p g V 


d Sv cop 

*37 - j FT3 Si + S ( z " z ) 


dz 


(0,z<) 



(L . z ' ) 


0 . 


(62) 


The boundary conditions on dgj/dz are a result of the conditions on V n . 
Implicit in (62) are conditions of continuity of g : at z = z' and discon- 
tinuity of amount j co/ p in the derivative of g : at z = z *. 

The result of applying these conditions to (62) is that 

co cos k ( z ' - L) cos k z 

g l < Z ’ 2 '> = 3 ¥ 13 sin k L ’ 2 '< z ' 

co cos k z 9 cos k ( z “ L) 
j T- k sinkL 7 ’ z > z' 


g v ( z > z ') 


- sin kz cos k ( z ' “ L) 
sin kL 


z < z 9 


- s in k ( z - L ) cos k z ' 
sin k L 


z > z' . 


(63) 
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The noise solutions thus become 


V„(«) = 


■-f 


sin k ( z - L) cos kz 1 
sin k L 


e( z ' ) dz' 


-r 

^ 7 


sin kz cos k (z' ~L) 

slrTkL 6(2 )d2 


I„(*) * i 






cos kz 9 cos k ( z ~ L) 
k sin k L 


e(z') dz' 


r 

* 7 


cos k ( z' - L) cos kz 
k sin k L 


e( z ' ) dz 




(64) 


The total output current (velocity) at z - L is obtained by adding I (L) and 


K (L) : 


" z 0 sinkL 

The spectral density S r (&i) of I(L) is obtained easily by using the fact that 
E s and e(z) are uncorrelated (i.e., E[E s e] = 0) and also using (33) and (36). 
Hence, upon multiplying (65) by its complex conjugate and taking the expected 


I. 


E + I cos kz ' e( z' ) dz # 


i 


(65) 
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value we obtain 



where is the signal spectral-density. 

Since the integrand appearing in (66) is non-negative it follows that in- 
creasing L, the length of the acoustic medium, results in a larger noise con- 
tribution to the output current (velocity) . This is to be expected because each 
elementary noise generator within the structure contributes noise independently 
to the output. Therefore, the longer the structure the more contributing noise 
sources present. 

In accordance with (26), the expected mean-square value of the output cur- 
rent is given by the integral (66) overall positive frequencies. If the signal applied 
at z = 0 is E s cos co s t then S s (w) = ttE s 2 s (w- <^ s ) and 




E l 

sin k L| 2 




N(«) da) , 


(67) 


where N(o)> is the second term in (66). By choosing L appropriately it is pos- 
sible to make the denominator of the first term in (67) small (if k(w s ) were 
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real, we could make the denominator vanish) and therefore amplify the signal, 
without amplifying all frequency components of the noise, simply by resonating 
the structure at the signal frequency. 

Signal Source and Matched Termination 

In this case we will assume that | a| « | . Under this condition Z Q , the 
characteristic impedance of the spin-phonon system, is approximately (/xp) 1/2 , 
which we take to be the terminating impedance. 

The boundary conditions are 

V(O) = E s 

V(L) * -I(L) • (/J.p) 1/2 , (68) 


which are also the boundary conditions satisfied by V 8 and I s . The boundary 
conditions satisfied by the noise are 

V n (0) = 0 

V n (L) = -^(L) • (.W>) l/a • (69) 

The Green’s function g x satisfies the boundary conditions 


£ii 

dz 


(O.z' ) 


0 , 


dgj 


[i 

°\ V2 1 


dz 

1 ( L, z ' ) 

= -j*(- 

ij g X 1 ( L,z ' ) ’ 

(70) 


39 



and the differential equations in (62). The solutions for g t and g v are 


g T (z- z') - 


-cos kz' e~^ kz 01 

k ‘ /X ’ 


z' < z 


e~i kz ' cos kz 

k /j. ’ 


z' > z 


g y (z, z') - coskz' e~^ kz , z' < z 


= j e •* kz ' sin kz , z ' > z , 


= ^KHa)) 


1/2 


(71) 


The total output voltage and current are given by 


V(L) = V s (L)+V n (L) = E s e-i kL + e“i kL f 

'o 


cos kz # e(z' ) dz ' 


I(L) = I s (L) + I n (L) 


( w >) 1/2 


e“ jkL 


e~ ikL 



cos kz' e(z' ) dz' . 


(72) 


In deriving (72) we have used the approximate value of Z 0 , stated earlier, and 
have set reflections equal to zero. 
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The spectral densities of (72) are 


S v O) - exp (2k : L) 


(f 


L \ 

S s (a>) + [ | | cos kz| 2 dz j y 


4Tlco 


m- 1 exp 



pco a 1 


1+a 


= exp (2k j L) S(co) + 


exp (4k r L) - 1 exp (2k , l) sin 2k R L 


4k, 


4k, 


Sj (co) 


S v (co) 

jj.p 



(73) 


kj, k R in (73) (the imaginary and real parts, respectively, of k) are given in 
terms of a by -b, a, respectively, of (54), As a check on (73) we note that if 
each T m = T > 0 then in the limit L- “, a - 0, in that order, we derive the results 
(56) for an infinite medium. 

For a signal spectrum concentrated at oj s 


E k 2 ] 


31 

2 


exp 


I CC 


N(o)) dco 


(74) 


where N(oi) is the second term in (73). 



If kj (co s ) is positive due to a population inversion of the energy- level pair 
whose frequency separation is co s , then we have amplification at the signal fre- 
quency (the phonon maser). A population inversion implies that for some m in 
(73), a.j andT ra are both negative. N(w) , however, remains positive, as it should. 

m 

It should also be noted in (73) that N(w) vanishes at those frequencies for which 
a vanishes. In our discussion of a. we pointed out that it was non- zero es- 

m 1 m 

sentially only in a bandwidth 2 /r^ m wide, centered at a% m . Hence, this is the 
band of frequencies which contribute to N(oj) . 
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